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ABSTRACT 

The  Nival  Postgraduate  School  computer  Simula t ion  model  QRAD  was 
used  to  study  sputtering  from  clean  and  nitrogen  reacted  tungsten,  and 
from  nitrogen  reacted  molybdenum,  for  (001)  surfaces  bombarded  by 
normal  Incidence  argon  Ions  In  the  energy  range  from  0.5  keV  to 
3.0  keV.  The  simulation  results  are  compared  to  experimental  data 
obtained  by  Winters,  who  concludes,  that  the  difference  in  the  cross 
sections  for  the  sputtering  of  nitrogen  from  W(001)  and  Mo(001)  is 
dominated  by  the  adsorbs te-substrate  mass  difference  since  the  other 
important  parameters  are  similar.  The  simulation  however  shows  little 
dependence  on  the  substrate  mass,  but  a  strong  dependence  of  the 
sputtering  cross  sections  on  the  distance  of  the  adatoms  above  the 
substrate.  The  nitrogen  position  is  required  to  be  higher  on  tungsten 
than  on  molybdenum. 
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I.  INTRODUCTION 


A.  OVERVIEW 

Sputtering  in  the  removal  of  surface  atoms  due  to  energetic 
per  tide  bombardment.  It  is  caused  by  collisions  between  incoming 
particles  and  atoms  in  the  selvage,  i.e.,  the  near  surface  of  solid 
material.  First  discovered  by  Grove  [1]  in  1853  and  Faraday  in  1854, 
it  took  about  50  years  until  the  physical  process  involved  as 
recognised  and  another  50  years  until  a  quantitative  theory  began  to 
develop.  Grove  [1]  noticed  the  disintegration  of  cathodes  in  glow- 
discharge  tubes  and  observed,  that  the  cathode  material  mas  deposited 
on  the  surrounding  glass  wells.  He  called  this  process  "cathode 
sputtering".  50  years  later,  Goldstein  l 21  presented  evidence  that  the 
sputtering  effect  mas  caused  by  positive  ions  of  the  discharge  hitting 
the  cathode. 

Stark  [3,4]  Introduced  a  so-called  'hot-spot'  model  of  sputtering 
(now  known  as  'evaporation  theory')  and  later  a  collision  theory,  which 
viewed  sputtering  as  a  sequence  of  binary  collisions  initiated  by  the 
bombarding  ion.  Stark  found  that  at  low  energies  the  yield,  Y, 
increases  with  increasing  ion  energy,  and  interpreted  this  result  on 
the  assumption  that  an  increasing  amount  of  energy  is  transferred  to 
the  terget  atoms  near  the  surface.  At  higher  energies,  increasing 
penetration  depth  of  the  ion  causes  a  relatively  smaller  sputtering 
yield  and  thus  a  plateau  la  the  yield  curve.  The  reason  for  this  is 
a  decrease  in  the  "slae”  of  the  ion  as  its  energy  increases.  Stark 
considered  his  hot-spot  model,  the  evaporation  of  target  material  from 


a  aicroscopica lly  small  region  with  high  local  temperature,  and  hi* 
collision  chaory  as  two  different  views  of  tha  sane  process. 

Von  Hippel  and  Blachschmidt  15],  who  developed  a  quantitative  for*  of 
Stark's  hot  spot  nodal,  aada  tha  first  attempt  to  formulate  a 
sputtering  theory  on  tha  basis  of  local  hast:  Subsequent  invest!** 

gators  found  that  this  model  is  inadquate.  Wehner  [6]  showed  that 
local  evaporation  alone  couldn' t  explain  sputtering  phenomena  from 
single  crystals.  In  1934  Lenar  and  Compton  [7]  suggested  that  binary 
collision  processes  dominate  in  light-ion  sputtering  and  local 
evaporation  is  dominant  in  hoevy-ion  sputtering.  Their  ideas  led  to 
the  'thermal  spike'  concept,  in  which  a  long  lived  high  temperature 
volume  presists  in  tha  target  after  the  collision  cascade  is  complete. 
Wehner' s  observations  revived  the  interest  in  the  collision  theory. 
Keywell  [8]  applied  neutron  transport  theory  to  Stark's  collision 
model.  Due  to  his  and  subsequent  calculations  by  Harrison  [9], 
probability  concepts  and  collision  cross  sections  were  introduced  as 
important  parts  of  the  sputtering  theory. 


B.  PHYSICAL  MEANING 

Sputtering  mechanisms  operate  in  the  near  surface  region  of  the 
target,  and  are  caused  by  a  collision  cascade.  The  processes  Involved 
are  similar  to  those  causing  radiation  damage  in  die  bulk  as  ter  la  1  of  a 
solid.  Sputtering  happens  far  from  tbs ram  1  equilibrium,  which  naans 
that  it  isn't  to  he  confused  with  evaporation.  Generally  an  incident 
particle  collides  with  atoms  of  tha  bulk  aaterlal  and  transfers  energy. 
If  the  trmasffered  energy  is  larger  than  dm  cohesive  energy  at  the 
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lattice  site  a  pr leery  recoil  atoa  la  created.  This  recoil  atom 


collides  with  other  target  a toss  and  distributes  the  energy  into  a 
collision  cascade.  The  condition  for  an  atom  to  be  ejected  is  that  the 
kinetic  energy  associated  with  the  normal  component  of  its  velocity  is 
larger  than  the  potential  energy  binding  it  to  the  surface.  The 
sputtering  yield  is  defined  as  the  average  number  of  target  atoms 
ejected  per  incident  particle. 

C.  LATER  WORK 

Although  considerable  discrepancies  exist  in  the  experimental 
results  produced  by  various  scientific  groups,  standards  have  developed 
which  define  the  basic  requirements  for  reproducible  and  reliable 
experiments.  This  development  is  shown  in  reviews  of  Behrisch  [10], 
Kaminsky  [111,  Carter  and  Colligen  [12],  and  Pleshivtsev  [13].  Over 
the  last  30  years  theoretical  efforts  have  been  increasingly  successful 
in  explaining  the  main  features  of  sputtering. 

The  analytic  theoretical  work  in  this  area  was  done  by  Sigmund 
[14].  He  presented  a  systematic  theory  of  sputtering  on  the  basis  of  a 
minimum  of  assumptions.  Basic  assumptions  concern  the  characteristics 
of  a  single  collision  (differential  cross  section,  inelastic  stopping), 
the  structure  of  the  medium  (random  or  crystalline),  the  structure  and 
binding  forces  of  the  surface,  and  some  factors  that  depend  on  specific 
experimental  arrangements  like  beam-target  geometry,  high  or  low  ion 
energy,  etc.  The  yield  calculation  consists  of  the  following  steps: 

1)  to  determine  the  amount  of  energy  deposited  by  energetic 
per  tides  (ion  and  recoil  atoms)  near  the  surface; 
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2)  to  convert  this  energy  into  a  number  of  low-energy  recoil  atoms; 

3)  to  determine  hew  many  of  these  atone  cone  to  the  surface; 

4)  to  select  those  atoms  that  have  enough  energy  to  overcome  the 
surface  binding  forces. 

Considering,  for  example,  the  backsputtering  yield  from  a  plane 
surface,  the  outward  current  of  the  target  atone  through  the  surface 
plane  at  x  »  0  is  found  to  be: 


J(E  ,&  )dE  d2a 
o  o  o  o 


r  dE  d2a 

fF  (e,  0,  o)  -H-5 -  Icos  el  — ° 

E  IdE  /dxl  4JI 

o  o 


where  E  is  the  energy  of  the  ion  impinging  at  an  angle  6  onto  a  plane 
surface,  Fq  is  the  deposited  energy  density,  T  is  the  number  of 
projectiles  per  unit  time,  0o  is  the  unit  solid  angle,  and  rm  is  a 
parameter  which  depends  on  H  and  the  scattering  cross  section.  The 
sputtering  yield  is  found  by  integrating  this  equation  over  E0  and 
0O  and  dividing  the  sputtered  current  by  the  incident  current  T. 

This  leads  to 

Y  -  A  Fd(E,  0,  0) 

with 

A  -  rm  f  dEo  t  d  (cos  e_)lcos  0.1  P  (E_.  0_) 

T  FW7ST  0  o  O  O 

o  o 


where  P(Eo,  ©„)  is  the  probability  that  an  atom  will  escape  from 

the  surfeca.  The  simplest  modal  for  the  surface  binding  of  a  metal  is 

based  on  a  planer  surface  barrier  U0.  In  this  case, 

1  for  Eq  cos^q  >  0o 
0  for  Bq  cos^o  <  U0 


P(*o»eo>  “ 
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which  leads  to  a  simple  form: 


4I2  N  C  U 

O  0 


with 

where  X0  2  24  and  a^  ■  0*219  A. 

D.  APPLICATION)  OF  SPUTTERING 

For  a  long  time  sputtering  was  regarded  as  an  undesired  effect 
which  destroyes  the  cathodes  and  grids  in  discharge  tubes.  In  fusion 
research,  plasma  contamination  by  metal  atoms  is  still  one  of  the  major 
problems.  Sputtering  also  may  cause  the  distructiou  of  diaphragms  and 
targets  in  accelerators  and  high-voltage  electron  microscopes. 

On  the  other  hand,  sputtering  has  become  an  indispensible  process 
in  modern  technology.  It  allows  the  removal  of  very  tightly  bound 
surface  layers  on  a  nearly  atomic  scale.  Sputtering  is  used  in 
sputter-ion  sources,  for  obtaining  reproducibly  clean  surfaces  and,  for 
example,  in  micro-machining  and  depth  profiling  of  thin  films. 

Sputtering  also  can  be  used  to  obtain  information  about 
processes  occuring  on  a  atomic  scale.  Many  scientists,  including 
Ringdon  and  Langmuir  (15),  Winters  (16),  and  Thglauer  and  Holland  (17), 
studied  the  sputtering  of  chemisorbed  gas  to  Investigate  mechanisms 
which  lead  to  the  emission  of  multi-component  materials.  This  thesis 
uses  a  computer  simulation  to  rnerlfy  one  assumption  made  by  Winters. 


E.  COMPUTER  SIMULATIONS 

Due  to  tneir  high  operating  speed,  computers  are  attractive  for 
calculations  concerning  sputtering  events.  A  computer  simulation  can 
never  replace  an  experiment,  but  by  matching  the  results,  a  simulation 
can  help  to  understand  and  identify  sputtering  mechanisms,  and  to 
interpret  experimental  results.  A  computer  simulation  is  not  a  theory, 
it  is  a  mathematical  tool  used  to  test  the  fundamentals  of  a  theory 
119]. 

Computer  simulations  can  be  divided  in  two  categories:  time— step 
and  event-store  models.  Event-store  programs  move  from  event  to  event 
skipping  the  intervening  time.  For  sputtering,  they  are  based  on 
binary  collision  theory.  These  fast,  complex  programs  are  useful  when 
the  physical  model  is  well  understood.  The  program  used  at  the  Ihval 
Postgraduate  School  is  a  time-step  model  based  on  simultaneous  multiple 
interactions  (MI).  MI  simulations  of  physical  systems  usually  require 
the  calculation  of  about  100  individual  ion  trajectories  (cascades), 
while  binary  collision  simulations  use  several  thousands.  This,  and 
the  fact  that  in  MI  simulations  mechanisms  are  easy  to  trace,  makes  the 
time-step  model  attractive. 

This  thesis  is  based  on  QRADS,  a  refined  version  of  QRAD,  the 
standard  program  in  use  at  the  Naval  Postgraduate  School  for  some  time, 
which  is  a  full-lattice  simulation  computer  program,  designed  to  model 
the  dissipation  of  an  incident  ion's  momentum  in  a  single  crystal 
target,  using  classical  mechanics.  The  program  has  the  ability  to 
generate  different  surfaces  of  several  crystal  structures.  To  obtain 
ra*cted  surfaces,  a  special  subroutine  places  adatoms  on  the  surface. 
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After  the  lattice  is  build,  the  potential  functions  are  calculated. 

The  iapotant  mathematical  manipulations,  concerning  the  trajectory,  are 
done  by  the  subroutine  INTEG  (integrator).  This  subroutine  also  tests 
for  the  termination  of  the  trajectory.  Finally,  subroutine  LEGIT 
determines  which  atoms  were  actually  ejected. 


II.  OBJECTIVES 


Some  experimental  results  exist  for  sputtering  from 
polycrystalline  BCC  structures  (191  but  only  a  few  for  clean  single 
crystal  surfaces  or  chemically  reacted  single  crystal  surfaces  [19]. 
Winters  [18]  investigated  the  sputtering  cross  sections  of  tungsten  and 
molybdenum  bombarded  by  argon  and  xenon  ions  in  the  energy  cange  from 
0.3  keV  up  to  5.0  keV.  Figure  1  shows  the  cross  sections  as  a  function 
of  incident  ion  energy.  Winters'  experiments  show  that  the  sputtering 
yield  of  nitrogen  tends  to  increase  as  the  atomic  weight  of  the 
surrounding  atoms  increases.  He  furthermore  concludes,  that  the 
difference  in  the  cross  sections  for  the  sputtering  on  nitrogen  from 
W(001)  and  Mo(001)  is  dominated  by  the  absorbs te-substrate  mass 
difference  since  the  other  important  parameters  are  similar. 

The  primary  objective  of  this  thesis  is  to  study  sputtering  from 


nitrogen  reacted  (001)  tungsten  and  molybdenum  surfaces  using  a 


computer  simulation  model  and  to  compare  the  results  with  the  exper¬ 
imental  data  obtained  by  Winters. 

The  incident  ion  will  be  Ar+  in  the  energy- range  from  0.5  keV  to 
3.0  keV. 

The  simulation  Includes: 

1)  Bombardment  of  a  clean  (001)  tungsten  surface  (used  for 
comparison). 

2)  Bombardment  of  a  nitrogen  reacted  (001)  tungsten  surface. 

3)  Bombardment  of  a  nitrogen  reacted  (001)  tungsten  surface  where 
the  atomic  mass  of  tungsten  is  replaced  by  the  atomic  mass  of 
molybdenum. 

4)  Bombardment  of  a  nitrogen  reacted  (001)  molybdenum  surface,  where 
the  nitrogen  is  closer  to  the  target  surface  plane  than  in  the 
cases  listed  above.  (The  nitrogen  is  assumed  to  be  in  a  4-fold 
symmetry  position  for  all  reacted  surfaces.) 


III.  THE  SIMULATION 


A.  THE  TARGET 

1.  Properties  of  Materials 

a)  Tungsten  (W,  Wolfram):  atomic  mass  183.85,  lattice  parameter 
3.165  A,  crystal  structure  BCC,  nearest-neighbor  distance 
2.74  A,  cohesive  energy  8.9  eV, 
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b)  Nitrogen  (N):  atomic  at as  14.0,  binding  energy  to  the  tungsten 
substrate  6.5  eV, 

c)  Molybdeaun  (Mo):  atomic  mass  94.94,  lattice  parameter  3.147  A, 
crystal  structure  BCC,  nearest-neighbour  distance  2.73  A, 
cohensive  energy  6.82  eV. 

2.  Location  of  the  Adatoms 

The  location  of  nitrogen  atoms  on  a  tungsten  or  molybdenum 
surface  is  not  known.  Clavenna  [21]  and  Adams  [22]  suggest  that  the 
nitrogen  is  single-bounded  in  a  4-fold  position  above  an  atom  in  the 
second  layer  of  the  tungsten  substrate.  (The  4-fold  position  is  also 
assumed  by  Winters.)  Figures  2a, b  show  side  and  top  view  of  the 
reacted  surface  where  the  radius  of  the  nitrogen  is  assumed  to  be  0.59 
A  [22].  (The  smaller  circles  are  the  nitrogen  atoms.)  Criffiths 
et«  al.  [23],  agree  with  this  position  at  a  fractional  coverage  (Of) 
of  0.5  but  mention,  that  at  this  coverage  the  nitrogen  is  unstable. 

They  recommend  study  at  the  most  stable  coverage  of  Of  *  0.4.  This 
requires  formation  of  Islands  of  nitrogen  atoms  on  the  surface,  (see 
Figs.  3a, b).  Harrison  et.  al.  [24],  propose  that  the  nitrogen  might  be 
below  the  first  layer  of  the  substrate  in  an  interstitial  site.  In 
this  case  some  second  layer  atoms  would  be  pushed  down  into  the 
substrate  to  set  equal  distances  to  all  nea rest-neighbors.  For  this 
simulation  the  nitrogen  is  assumed  to  be  in  the  4-fold  position  for 
both  substrate  materials.  For  the  nitrogen  reacted  tungsten  surface 
the  location  above  the  substrate  is  chosen  so  that  the  distance  to  all 
five  nea  rest-neighbors  is  equal.  The  distance  above  the  surface  is 
0.504  LU  (where  1.0  UJ  ■  1. 58  A  ■  0.5  lattice  parameters,  see 


Fig.  4) 


(Lattice  units  were  introduced  in  the  program  so  that  rescaling  to  a 
new  target  material  can  be  dona  without  major  reprogramming  and  so  that 
results  from  different  substrates  can  be  compared  directly.)  Runs  also 
were  made  with  the  tungsten  mass  replaced  by  the  molybdenum  mass  to 
examine  the  ass  effect  proposed  by  Winters  118],  discussed  in  chapter 
2.  On  the  nitrogen  reacted  molybdenum  surface  the  nitrogen  Is  placed 
closer  to  the  surface  layer  (see  Fig.  4).  The  distance  above  the 
surface  is  0.243  LU  (where  1.0  LU  »  1.57  A  •  0.5  lattice  paxametrs). 

3.  The  Target  Size 

The  inability  of  the  lattice  of  target  atoms  to  contain  the 
entire  cascade  generated  by  the  incident  ion  is  called  '  a  failure  in 
containment1 •  The  size  of  the  target  determines  the  containment  of  the 
trajectory.  If  an  Increase  in  target  size  doesn't  change  the  results 
of  the  computations,  the  trajectory  is  effectively  contained  for  the 
purposes  of  the  computations.  (Unnecessarily  large  targets  waste 
computer  resources. ) 

The  minimum  target  size  for  the  ion  energy  of  0. 5  keV  was  found 
to  be  19x8x19,  where  the  numbers  represent  the  numbers  of  planes  In  the 
x,y  and  z  directions.  40  nitrogen  atoms  (equivalent  to  a  fractional 
converage  of  ■  0.4)  are  placed  as  described  above  (see  Fig.  S; 
solid  lines  mark  the  smaller  target).  The  minimum  target  size  for  ion 
energies  between  1.0  keV  and  3.0  keV  was  found  to  be  23x8x23  with  a 
fractional  coverage  of  9f  ■  0.416,  equivalent  to  60  adatoms.  The 
difference  in  coverage  does  not  significantly  affect  tbs  result. 

Winters  [18] ,  determined  the  coverage  in  his  experiments  by  Auger 
spectroscopy,  and  found  that  the  Auger  intensity  is  a  nonlinear 
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,2 


Thus  all 


function  of  coverage  above  0  •  6  *  101**  atone  per  cn  . 
his  measurements  vers  taken  below  this  value* 

The  coverage  chosen  for  the  s inula t ion  is  equivalent  to  4. 9  x  10^ 
atoms  per  cm2. 


B.  POTENTIAL  FUNCTIONS 
1*  In  General 

A  primary  uncertainty  in  sputtering  models  and  theory  is  the 
form  of  the  interatomic  potential  functions.  No  single  potential 
function  has  yet  been  deduced  which  fits  every  sputtering  experiment. 
The  functions  used  for  this  simulation  are  Morse,  and  Moliere  [25] 
joined  with  a  cubic  spline  [24].  The  Morse  potential  function  is  of 
the  form: 

V(r)  -  De  exp[-2a(r-r0)]  -  2De  exp[-a(r-r0)] 
where  De  is  the  well  depth,  a  is  a  scale  factor  of  the  potential 
function,  and  rQ  is  the  equilibrium  separation.  For  r  >  rQ  the 
Morse  potential  is  attractive  and  for  r  <  rQ  repulsive. 

The  Moliere  potential  function  is  of  the  form: 

V(r)  -  (ZiZ2e2/r)(0.35  exp[-0.3r/a]  +  0.55  exp[-1.2r/a] 

40.10  exp[-6.0r/a]) 

where  'a*  is  the  Flrsov  screening  length.  The  Moliere  function  is  a 
pure  repulsive  potential  function. 

The  cubic  spline  is  of  the  form: 

V(r)  "  %  ♦  Cj^  r  ♦  Cjr2  Cjr3 

where  the  coefficients  are  determined  by  tbs  program.  The  cubic  spline 
is  necessary  to  combine  Morse  and  Moliere  functions  into  a  single 
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composite  potential.  For  hard  collisions,  at  close  approach  in  the 
keV-range,  the  Morse  potential  bends  towards  the  Y-axis  and  thus  cannot 
be  physically  correct,  because  all  potential  functions  oust  approach  “ 
as  r  ♦  0  like  1/r.  The  cubic  spline  is  inserted  where  the  Morse  and 
the  Moliere  functions  overlap,  thus  the  contained  function  is  a  Moliere 
to  the  left  and  a  Morse  to  the  right  of  the  spline.  In  general  the 
spline  shouldn't  be  longer  than  0.05  LU  to  get  smooth  potential  and 
force  curves. 

2.  Crystal  Potential  Functions 

The  composite  potential  functions  are  used  for  the  bulk 
materials  as  well  as  for  the  reacted  surfaces.  The  Moliere  function  is 
a  standard  Moliere  in  each  case.  The  Morse  parameters  have  physical 
meaning:  a,  used  as  an  adjustable  parameter,  is  chosen  so  that  the 
slope  of  the  Morse  matches  that  of  the  Moliere  function.  For  the 
substrate  materials  De  is  chosen  to  generate  individual  atom  energies 
in  the  center  of  the  target  equal  to  the  cohesive  energy,  i.e.,  8.9  eV 
for  tungsten  and  6.82  eV  for  molybdenum.  For  the  reacted  surfaces  Dg 
has  to  match  the  binding  energies  between  the  nitrogen  and  the  substate 
i.e.,  6. 5  eV  in  each  case.  Because  the  potential  function  is  truncated 
after  one  nearest-neighbor  distance,  r0  is  greater  than  the 
separation  between  atone  in  the  substrate  materials.  For  the  reacted 
surfaces,  rQ  is  the  nearest-neighbor  distance. 

The  cubic  spline  la  bounded  by  the  values  ft&  and  ftg.  The 


table  below  shows  the  parameters  for  the  crystals  and  adatoms 


Table  1:  Parameters  for  Crystal  Potential  Function. 


«(A'1) 

r«(A) 

De(eV) 

RX(A) 

RgCA) 

see  Figs. 

w  s* 

1.20 

2.849 

1.335 

1.83 

1.90 

6,10 

W  -N 

1.65 

2.380 

1.300 

1.14 

1.14 

7,10 

Mo-Mo 

1.50 

2.800 

0.997 

1.24 

1.30 

8,10 

Mo-N 

2.20 

1.950 

1.612 

0.91 

0.97 

9.10 

Figures  6,  7,  8, 

9  show 

the  Morse  and 

Mollere 

functions  and  Figure 

10  shows  the  combined  functions  generated  by  the  program. 

3.  Gas  Phase  Substrate  Atom  Potential  Functions 

The  gas  phase  potential  functions  are  pure  Morse  functions 
because  no  hard  collisions  are  expected.  The  parameters  are  listed 
below: 


Table  2:  Parameters  for  Gas  Phase  Substrate  Atom  Potential  Functions. 


o(A“1) 

r«(A) 

D.(eV) 

see  Figs. 

W  -W 

1.20 

2. 74 

8.900 

11 

W  -N 

1.65 

1.96 

6.721 

11 

Mo-Mo 

1.50 

2.73 

6.820 

12 

Mo-N 

2.20 

1.96 

5.100 

12 

4.  Nitrogen-Nitrogen  Potential  Functions 

The  N-N  potential  function  In  the  gas  phase  as  well  as  In  the 
crystal  phase  la  a  pure  Morse  functions  tilth: 

a  -  2.7  A’1 

r«  •  1.098  A 

D9  -  7.373  eV  (crystal  phase) 

D*  -  9.834  eV  (gas  phase) 

see  Figures  11,  12  for  the  gas  phase  potential  functions. 

3.  Argon-Tungsten/ Molybdenum  and  Arson- Nitrogen  Potential 
Functions 

Bscs use  it's  not  likely  that  argon  Ions  reset  with  nitrogen 
tungsten  or  molybdenia,  these  potential  functions  here  the  standard 
Hollers  form,  (see  Fig.  13). 
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C.  THE  IMPACT  AREA 


The  Impact  area  for  each  simulation  is  a  rectangle  containing  a 
grid  of  104  impact  points,  (sea  Fig.  14),  where  the  length  of  the  sides 
is  3.165A  “  2.0  LU  for  tungsten  and  3.147  A  -  2.0  LU  for  molybdenum. 

The  origin  of  the  rectangle  is  placed  at  the  following  locations: 


Table  3:  Origns  of  the  Impact  Area 


target  sise 


locations 


23x8x23 

19x8x19 

23x8x23 

19x8x19 


N-raacted  (10,10)  (10,12)  Fig.  5 
N-reacted  (08,08)  (08,06)  Fig.  5 
dean  (11,11)  Fig.  5 
dean  (09,09)  Fig.  5 


On  the  NLtrogen-raected  surface,  two  impact  point  sets  are  neces¬ 
sary  to  assure  complete  converage. 


IV.  DISCUSSION  AM)  RESULTS 


All  output  hus  analysed  by  a  separate  analyser  program,  AN83,  which 
provides  a  nuawrical  analysis,  and  by  the  program  PLOTAN,  providing  a 
graphical  analysis.  The  latter  is  discussed  first. 

A.  GRAPHICAL  ANALYSIS 
1.  Spot  Pattern 

The  spot  pattern  plots  show  the  polar  and  aslmuthal  angular 
distribution  of  the  ejected  atoms,  (see  Figs.  13,  16,  17).  The 
distances  of  the  crosses  form  the  center  is  a  aaasure  of  the  polar 
angle  of  the  ejected  atoms,  l.e.,  a  eross  far  out  corresponds  to  a 
large  ejection  angle  measured  from  the  surface  normal.  Tbs  spot 
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pattern  obtained  for  the  substrate  materials  (W,  Mo,  Mo-mss)  show  a 
four-fold  systry,  where  'Mo-mass'  means,  that  in  the  program  the 
tungsten  sees  ms  replaced  by  the  molybdenum  mass  and  all  other 
parameters  were  unchanged. 

The  spots  at  45® ,  135®,  255*  and  315*  point  in  the  direction  of 
the  naa rest-neighbors.  Figure  15  shows  the  spot  pattern  for  molybdenum 
atoms  ejected  from  the  substrate  after  bombarding  a  nitrogen-reacted 
molybdenum  surface  by  2.0  keV  argon  ions.  This  figure  is  typical  of 
all  spot  pattern  obtained  from  substrate  materials. 

The  spot  pattern  obtained  from  the  adatoms  show  no  symmetry. 
Figures  16,  17  are  representative  for  ejected  nitrogen  atoms.  Figure 
16,  the  case  where  the  nitrogen  is  placed  further  above  the  surface, 
shows  a  heavier  concentration  of  ejected  atoms  at  low  angles  than 
Figure  17. 

2.  Energy  Distribution 

Figure  18  shows  typical  energy  distributions  for  atoms 
sputtered  from  the  substrate  in  the  ion  energy  range  from  0.5  keV  to 
3.0  keV.  The  envelopes  peak  near  7.0  eV. 

The  energy  distribution  of  the  sputtered  nitrogen  atoms  depends 
on  the  location  of  the  adatom.  If  the  nitrogen  is  close  to  the  surface 
layer,  the  distribution  peaks  at  about  2.5  eV/atom,  (see  Fig.  19).  If 
the  nitrogen  is  further  away,  as  in  the  case  of  the  tungsten  or 
molybdenum-mss  substrate,  the  distribution  is  much  broader,  and  the 
maximum  is  roughly  15  eV/atom,  (see  Fig.  20). 
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3.  Energy  Distribution  >100  ggec 

The  energy  distribution  for  particles  ejected  from  the 
substrate  after  100  fsec  (1.0  fsec  "  10“^  sec)  is  similar  for  all 
materials,  but  with  the  difference  that  the  fraction  for  molybdenum  or 
molybdenum-mass  is  less  than  for  tungsten,  (see  Figs.  21,  22,  23). 

More  low  energy  atoms  are  ejected  after  100  fsec  in  the  case  of  the 
heavier  substrate  material.  This  could  indicate  a  mass  effect. 

Energy  distributions  for  the  ejected  nitrogen  atoms  are  shown 
in  Figures  24  and  2S.  If  the  nitrogen  is  in  the  lower  position  there 
were  no  atoms  ejected  with  energy  greater  than  7.0  eV.  This  behavior, 
is  seen  throughout  the  whole  energy  range,  (0.5  kaV  to  3.0  keV). 
Although  the  total  number  of  atoms  is  low,  (see  Fig.  25),  and  therefore 
the  statistics  are  poor,  the  difference  between  Fig.  24  and  Fig.  25  may 
be  indicating  a  distance  effect. 

4.  Energy  Versus  Time. 

Figures  26,  27,  28  and  29  show  that  no  high  energy  particles 
are  ejected  later  than  100  fsec.  The  spectrum  for  substrate  materials 
(Fig.  26)  is  wider  than  that  for  the  ejected  nitrogen  atoms. 

Figure  28  shows  that  hardly  any  nitrogen  atoms  are  ejected 
after  100  fsec.  The  difference  between  Fig.  27  and  29  and  Fig.  28  on 
the  other  hand  shows  again,  that  there  may  be  a  distance  effect. 

5.  Ejection  Time  Distribution 

The  ejection  time  distribution  is  essentially  the  same  for  all 
substrate  materials.  For  high  energies  the  distribution  is  lower  and 
broader,  (see  Figs.  30  aad31).  The  distribution  peaks  at  about  80 
fsec. 
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6.  Atoms  £gr  Single  Ion 

The  atoms  per  single  ion  diagrams  are  similar  throughout  all 
runs.  The  spectrum  broadens  with  increasing  energy,  (see  Figs.  38,  39, 
40  and  41). 

7.  Yield  per  Single  Impact  Point 

The  yield  per  single  impact  point  diagrams  are  also  similar  for 
all  runs,  (see  Figs.  42  and  43),  where  the  open  circles  represent  no 
yield.  The  diagrams  are  combined  for  the  locations  (10,  10)  and  (10, 
12),  i.e.,  the  upper  left  corner  of  each  rectangular  box  represents  the 
location  (10,  10).  Comparison  of  the  centers  of  the  combined  squares 
shows,  that  whenever  the  ions  hit  a  substrate  atom  the  yields  are 
larger.  The  center  of  each  right  square  is  occupied  by  a  nitrogen 
atom.  The  nitrogen  acts  like  a  'buffer'.  The  open  circles  are  caused 
by  channeling.  Figures  42  and  43  show  that  with  decreasing  ion  energy 
the  areas  where  -channeling  occurs  decrease.  The  reason  for  this  is 
that  for  decreasing  energies  the  cross  sections  increase  because  the 
impact  energy  appears  in  the  denominator. 

8.  Yield  per  Surface  Atom 

Figures  44,  45  show  the  yield  obtained  from  the  surface  layers 
for  both  target  sites.  Almost  no  atoms  are  ejected  from  the  edges  of 
the  crystal.  This  indicates  that  the  crystals  are  large  enough  to 
contain  the  trajectories. 

B.  NUMERICAL  ANALYSIS 

The  numerical  analysis  done  by  the  analyser  program  AN  83  provides 
the  following  nitrogen  yields,  (see  also  Fig.  46): 
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Table  4:  Nitrogen  Yields  from  Simulations. 


substrate: 

W(001) 

Mo-mass(OOl) 

Mo(001) 

energy 
in  keV 

yields 

yields 

yields 

0.5 

2.36 

2.04 

0.60 

1.0 

2. 25 

2.29 

0. 68 

2.0 

1.80 

2.19 

0.51 

3.0 

1.37 

1.87 

0.51 

Experimental  yields  in  the  ion  energy  range  from  0.5  keV  to  3.0  keV 
are  known  only  for  polycrystalline  tungsten  and  molybdenum  surfaces 
[26].  Figure  47  shows  the  experimental  metal  yields  obtained  from 
polycrystalline  tungsten  and  molybdenum  as  well  as  the  yield  from  clean 
U(001)  obtained  from  the  simulation.  The  yields  are  listed  below: 

Table  5:  Yields  from  Clean  Metal  Surfaces. 

Experiment  Simulation 


energy 
in  keV 

W(poly) 

Mo(poly) 

W(001) 

0.5 

0.58 

0.57 

2.31 

1.0 

0.88 

0.80 

2.94 

2.0 

1.30 

1.10 

2.40 

3.0 

1.50 

1.30 

2.18 

Some  data  are  known  for  5.0  keV  Ar+  bombardment  of  molybdenum 
[20]  at  normal  Incidence  onto  the  (001)  planes  and  onto  polycrystalline 
molybdenum.  The  yield  from  polycrystalline  Mo  is  found  to  be  larger 
than  from  Mo(001).  Thus  Roosendaal  [20]  concludes  that  for  a  BCC 
crystal  the  same  result  is  found  as  for  FCC  structures:  the  sputter¬ 
ing  yield  for  less  closed  packed  planes  is  less  than  the  yield  from 
polycrystalline  targets,  i.e.,  Y(Foly)  >  Y(001).  The  results  of  this 
simulation  show  a  yield  of  about  twice  the  yield  from  polycrystalline 
tungsten.  Winters  [16]  and  Vfehner  [27]  who  bombarded  W(001)  surfaces 
in  the  low  energy  zange  with  various  ion  types  obtained  the  results 
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shown  below: 


Table  6:  Experimental  Yields  from  Clean  Metal  Surface  for  Ion 

Energy  of  0.5  keV. 


Neon 

Argon 

Krypton 

Xenon 

Winters : 

0.25 

-  .1 — 

0.91 

1.21 

Wehner : 

0.27 

— 

0.92 

1.02 

The  yield  Is  Increasing  with  Increasing  Ion  mass  and  should  be 
about  0.6  for  0.5  keV  Ar+  Ions.  The  result  of  this  simulation  Is 
higher  than  this  value  by  a  factor  of  about  4.  This  shows  that  the 
(Ar-W)  potential  (a  standard  Mollere)  Is  too  hard  In  this  energy 
range. 

Winters  [18]  shows  his  results  as  cross  section  versus  Ion  energy 
(see  Fig.  1).  The  nitrogen  cross  sections  are: 

Table  7:  Experimental  Nitrogen  Cross  Sections. 

cross  sections  x  10”^,  (cm^) 


energy 
in  kev 

substrate 

W(001) 

Mo(001) 

0.5 

0.91 

0.65 

1.0 

0.98 

0.72 

2.0 

1.00 

0.75 

3.0 

0.94 

0.73 

Winters  assumes  that  the  nitrogen  sputtering  rate  is  given  by 
"rN  “  dOjj/dt  -  OnOnV+, 

where  0))  Is  the  nitrogen  atom  concentration  (atoms/ cm^),  oN 
Is  the  sputtering  cross  section  (ca^),  v+  the  Ion  flux  (ions/cm^- 
s),  and  t  the  time  In  seconds.  The  product  tv+  is  the  Ion  dose 
(fluence)  While  the  product  gives  the  sputtering  yield  in 
nitrogen  atoms/ ion. 


Using  this  relation  (Y  •  QjfJ }j)  the  cross  section  obtained  from 
the  simulation  are: 


Table  8: 

Simulation 

Nitrogen  Cross 

Sections. 

cross 

sections  x 

10"15,  (cm2) 

energy 
in  keV 

substrate: 

W(001) 

Mo-mass(OOl) 

Mo(001) 

0.5 

4.77 

*.12 

1.20 

1.0 

4.55 

1.36 

2.0 

3.  64 

*.  *»2 

1.02 

3.0 

2.77 

1.  78 

1.02 

Figures  48  and  49  show  Winters  cross  sections  In  coape r Ison  to  the 
cross  sections  obtained  from  tbs  simulation.  In  the  case  where  the 
nitrogen  is  in  the  lower  position  (Mo)  (Fig.  48),  the  simulation  cross 
sections  are  in  a  better  agreeaent  with  the  experimental  results.  But 
even  In  this  case  the  simulated  cross  sections  are  still  too  large. 

This  suggests  that  the  ( Ar-N)  potential  (also  a  standard  Moliere)  is 
too  hard.  Winters'  experiments  showed  that  the  sputtering  yield  of 
nitrogen  tends  to  Increases  as  the  atomic  weight  of  the  substrate 
increases.  He  suggests  that  the  difference  in  the  cross  section  for 
the  sputtering  of  nitrogen  from  W(001)  and  Mo(001)  is  dominated  by  the 
adsorbate-substrate  mass  difference  since  the  other  Important 
parameters  are  similar.  The  simulation  shows  that  the  distance  between 
the  nitrogen  and  the  substrate  is  more  important  than  the  mass 
difference. 

Winters  furthermore  states  that  the  shape  of  the  yield  curves  for 
the  chemisorbed  gas  and  the  substrate  material  is  different  because  the 
ratio  o(gas)/o(substrate)  decreases  with  increasing  ion  energy. 

Winters  used  o (substrate)  cross  sections  obtained  from  polycrystalline 
W  and  Mo  asstalng  a  surface  coverage  of  0  •  10+^  atoms/ cm2.  The 
following  ratios  were  calculated  from  the  simulation:  oN(W)/oW(clean) , 


o (Mo-mass) /oW(clean) ,  oN(Mo)/aW(clean) .  Using  aW(clean)  in  the 
latter  case  seems  reasonable,  because  comparision  of  yields  from 
polycrystal-  line  W  and  Mo  as  well  as  the  yields  obtained  from  the 
substrates  are  not  too  different,  i.e.,  the  yields  from  Mo(clean) 
shouldn’t  divlate  much  from  W(dean).  ©  for  W(001)  is  9.98  x  10+^ 
atoms/cm^.  The  ratios  are  shown  below  (see  also  Fig.  49). 

Table  9:  Cross  Section  Ratios. 


Experiment 

Simula  tion 

oN(W) 

oN(Mo) 

oN(W) 

a N(Mo-ma88) 

a  N(Mo) 

energy 
in  keV 

aW(poly) 

aMo(poly) 

o w( clean) 

aW(clean) 

oW( clean) 

0.5 

1.57 

1.14 

2.06 

1.78 

0.52 

1.0 

1.11 

0. 90 

1.54 

1. 57 

0.46 

2.0 

0.77 

0.68 

1.52 

1.84 

0.43 

3.0 

0.63 

0.56 

1.27 

1.73 

0.47 

In  comparison  to  Winters'  results.  Figure  30  shows  the  ratios  for 
nitrogen  in  the  lower  position  and  Figure  51  for  nitrogen  in  the  higher 
position.  (Mo-mass  not  shown.)  The  values  obtained  from  the  simula¬ 
tion  show  the  general  trend  observed  by  Winters  for  the  W  and  Mo 
substrate,  but  Indicate  that  the  location  of  the  high  nitrogen  is 
perferable  for  the  tungsten  substrate  and  low  for  the  molybdenum 
substrate.  Even  with  these  choices,  the  difference  between  the 
nitrogen  locations  is  too  large.  This  is  also  seen  by  comparing 
oN(W)/oN(Mo)  between  Winters  data  and  the  data  obtained  from  the 
simulation.  The  ratios  are  shown  below,  (see  also  Fig.  32): 

Table  10:  Cross  Section  Ratios. 

energy  Experiment  Simulation 

in  keV  <JN(W)/cN(Mo)  oN(W)/oN(Mo-t»ss)  oN(W)/oN(Mo) 


0.5 

1.40 

1.16 

3.98 

1.0 

1.36 

0.98 

3.35 

2.0 

1.33 

0.82 

3.57 

3.0 

1.29 

0.73 

2.72 
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As  mentioned  in  section  A7,  with  decreasing  ion  energy  the  cross 
sections  increase  and  the  areas  where  channeling  occur  decrease.  This 
is  also  seen  comparing  the  reflected  argon  yields,  listed  below: 


table  11:  Simulation  Argon  Yields. 


clean 


energy 
in  keV 

substrate  W(001) 

0.5 

0.40 

1.0 

0.27 

2.0 

0.12 

3.0 

0.07 

nitrogen 

reacted 

Mo-mass(OOl) 

W(001) 

Mo(001) 

0.30 

0.43 

0.24 

0.11 

0.22 

0.13 

0.07 

0.09 

0.05 

0.07 

0.05 

0.03 

Due  to  the  increase  in  cross  sections,  more  ions  are  reflected  at 
low  energies.  Harrison  (19)  says  that  about  85%  of  the  ejected  atoms 
(bulk  material)  from  clean  targets  come  from  the  first  layer.  This 
simulation  shows  that  at  least  802  of  the  ejected  atoms  come  from  the 
first  layer.  The  yield  from  deeper  layers  varies  between  192  and  52, 
where  the  yield  from  below  the  second  layer  is  almost  neglible.  The 
yields  from  below  the  first  one  are  shown  below: 


Table  12:  Yields  from  Second  or  Deeper  Layers  (Simulation), 
clean  nitrogen  reacted 


energy 
in  keV 

substrate 

W(001) 

Mo-mass(OOl) 

W(001) 

Mo(001) 

0.5 

112 

122 

162 

52 

1.0 

102 

142 

172 

82 

2.0 

162 

162 

192 

122 

3.0 

172 

182 

172 

162 

Some  of 

the  ejected  atoms,  which 

have  similar 

energy, 

velocity  and 

position,  combine  above  the  surface,  forming  multimers.  Moat  of  the 
multimers  observed  are  dimers:  (W-W,  (W-N),  Mo-Mo),  (Mo-N).  No  (N-N) 
dimers  were  found.  The  numbers  of  dimers  observed  are  listed  below: 


2 


table  13:  Dimers  from  Clean  Tungsten  (Slmultlon). 


energy 

W  - 

in  keV 

0.5 

40 

1.0 

53 

2.0 

18 

3.0 

22 

Table  14: 

Dimers  from 

Nitrogen  Reacted 

Surfaces 

(  Simula  tion) 

energy 

W  -  N 

Mo  -  N 

W  -  W 

Mo  -  Mo 

in  keV 

0.5 

50 

9 

18 

18 

1.0 

65 

12 

23 

31 

2.0 

27 

15 

24 

15 

3.0 

14 

8 

19 

10 

The  table  below  shows  the  number  of  linked  atoms  (taken  from  all 
mil timers) . 


Table  15:  Linked  Atoms 

(Simulations). 

energy  substrate 

U 

Mo-tm 

ss 

Mo 

in  keV 

(WW) 

(WN) 

(MoMo) 

(MoN) 

(MoMo) 

(MoN) 

0.5 

33 

71 

56 

85 

33 

18 

1.0 

41 

82 

74 

92 

54 

29 

2.0 

72 

66 

37 

64 

27 

18 

3.0 

42 

37 

42 

65 

29 

19 

In  the  case  of  the  molybdenum  substrate  the  cumber  of  Mo-N  links 
exceeds  the  number  of  Ho-Mo  links  throughout  the  whole  energy  range. 
In  the  case  of  the  tungsten  substrate  the  nueber  of  linked  W-W  is 
higher  for  3  keV  end  2  keV. 
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V.  SUMMARY  Ate  RECOMMENDATIONS 


The  results  of  the  simulation  show  that  the  lon-substrate 
potential  functions  are  chosen  too  hard*  Further  work  should  be  done 
to  reduce  the  yields  for  the  substrate  materials.  The  graphical 
analysis  shows  that  there  may  be  mass  as  well  as  distance  effects,  but 
the  results  of  the  simulation  clearly  show  that  the  adsorbate— substrate 
distance  is  more  important  than  the  adsorbate— substrate  mass 
difference,  when  the  nitrogen  sits  in  the  higher  position  on  the 
substrate  with  tte  higher  atomic  mess.  After  adjusting  the  substrate 
yields,  futher  work  should  be  done  to  find  the  right  distance  between 
the  adatoms  and  the  substrate. 
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Fig.  1.  Sputtering  Cross  Sections. 
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a)  Plan  via?;,  illustrating  domain  and 
boundary  structure. 

b)  Cross  section  through  the  domain 
along  AAl 

Island  formation  of  a  nitrogen  reacted 
tungsten  surface. 


b)  Topvlew 

A  x 

first  layer 
second  layer 

'A' 

Small  circles  represent  nitrogen  atoms. 

" 

Fig.  4*  Location  of  nitrogen  atoms  on  tungsten 
and  molybdenum  substrates. 
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Small  circles  represent  nitrogen  atoms 
Solid  lines  mark  smaller  target. 
Numbers  represent  layers. 

Fig*  5*  Topview  of, reacted  surface. 
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potential  energy,  ( eV ) 


Separation,  (A) 

Fig.  8.  Interatomic  potential  functions. 
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potential 


Fig.  9.  Interatomic  potential  functions. 
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otential  ener 


Fig.  10.  Interatomic  potential  functions. 
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eparation,  (LU) 


SPOT  PATTERN 

(obtained  for  Mo-substrate) 


3.0  KEV  W1001)  N* 

SPOT  PATTERN 
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3.0  KEV  M0(001)  *N* 


SPOT  PATTERN 


ENERGY  DIS I RIBU1  ION  MOIOO  I)  MO' 
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ENERGY  (EV)  ENERGY  (EV) 

Energy  distribution  for  ejected  Mo  atoms  from  a  reacted  surface 


ENERGY  DISTRIBUTION 


NOllOVfel 


0.0  0.0  to  0  IC  O  20.0  2C.0 

ENERGY  (EV) 

Fig.  20.  Energy  distribution  for  ejected  N  atoms  from  reacted  surfaces. 
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bution  for  atoma  (substrate),  ejected  after  100  fsec 


M.  (100)M 


0.6  KEV  MCM001)  *N* 


TIME  DISTRIBUTION 


7?.  Election  time  distribution  for  N 


Fig.  33.  Ejection  time  distribution  i'or  H  atoms 


distribution  for  !J  atoms 


stribution  for  N  atom 


FRACTION 


O.B  KEY  WtOOl)  *N* 


ATOMS  PER  SINGLE  ION 


0  10  20 

ATOMS  PER  SINGLE  ION 

Fig.  33.  ASI  for  ejcctod  N  atoms. 
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to  KEV  wiooo  *N* 


ATOMS  PER  SINGLE  ION 


O  10  so 


ATOMS  PER  SINGLE  ION 

Fig.  39*  ASI  for  ejected  X  atoms. 
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FRACTION 


2.0  KEV  W(001)‘N* 


..  ATOMS  PER  SINGLE  ION 


— i  -  i 

O  10  M 

ATOMS  PER  SINGLE  ION 

Fig.  VO.  AS I  for  ejected  N  atoms. 
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10  so 

ATOMS  PER  SINGLE  ION 

V.  atoms. 
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Tig,  44.  Held  per  surface  atom. 
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Fig.  45«  Yield  per  surface  atom. 
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Fic.  47.  Simulation  yields  from  V.'(OOI)  and  experimental  yields 
from  polycrystalline  'I!  and  Mo. 
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